Resonant control of polar molecules in an optical lattice 
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We study the resonant control of two nonreactive polar molecules in an optical lattice site, focussing on the 
example of RbCs. Collisional control can be achieved by tuning bound states of the intermolecular dipolar 
potential, by varying the applied electric field or trap frequency. We consider a wide range of electric fields 
and trapping geometries, showing that a three-dimensional optical lattice allows for significantly wider avoided 
crossings than free space or quasi-two dimensional geometries. Furthermore, we find that dipolar confinement 
induced resonances can be created with reasonable trapping frequencies and electric fields, and have widths that 
will enable useful control in forthcoming experiments. 

PACS numbers: 03.65.Nk, 34.10,+x, 34.50.Cx 



Ultracold gases of polar molecules are of interest for their 
long-range dipolar interactions, which give them unique ap- 
plications in areas such as many-body phases [1], quantum in- 
formation [2], and precision measurement [3-5]. Cold gases 
of LiCs [6] and RbCs [7] have been formed with tempera- 
tures T < 1 mK, and work continues to produce degenerate 
gases [8, 9]. Since the creation of a near-degenerate gas of 
40 K 87 Rb with T < 1 jiK [10], a number of studies have been 
done of its collision properties [11-14]. KRb has an exother- 
mic reaction producing K 2 + Rb 2 , which occurs with almost 
unit probability when two molecules are sufficiently close. 
This allows a simple description of the collision properties in 
terms of universal physics [15-22]. Such techniques should 
also apply to other species with reactive collisions, and to the 
quenching of any vibrationally excited molecule. For many 
studies it is therefore desirable to keep molecules separated, 
for example by confining the gas in a three-dimensional ( '3D' ) 
lattice [23] or in a quasi-two dimensional ('2D') geometry 
with the molecules polarized perpendicular to the plane ( 'side- 
by-side') [14-17, 24, 25]. This reduces the likelihood of 
molecules approaching each other along the attractive ' head- 
to-tail ' path which is available in 3D. 

In contrast to reactive molecules such as KRb, ground state 
NaK, NaRb, NaCs, KCs and RbCs are not reactive, and so are 
available for experiments on longer timescales and at higher 
densities where control of elastic collisions is useful. The long 
range dipole-dipole interaction between two molecules pro- 
duces an anisotropic potential which is capable of supporting 
bound states [26] . Tuning these bound states around a colli- 
sion threshold with an electric field allows resonant control of 
the interactions [18, 27, 28], in analogy to the magnetic and 
optical control that has been so useful for neutral atoms [29] . 
Because three-body recombination can still occur [30], iso- 
lating a pair of molecules in an optical lattice site provides 
an ideal, loss-free environment for studying the two-body en- 
ergy spectrum. Such a scenario is analogous to several ex- 
periments performed on atom pairs [31-34]. Optical lattices 
have also been used to tune atomic collisions through con- 
finement induced resonances [35-39] (CIRs), which depend 
on the scattering length being comparable to the characteris- 



Quantity 


Definition 


Value (fl ) 


Mean scattering length 


* = WiW^ 2m ^ in2)llA 


233.5 


Confinement length 


4o = ^h/(2m r aj) 




cj/ln = 1 kHz 




5728 


o)/2tt = 50 kHz 




810.5 


Dipole length 


a^ - m^ji 1 /ft 2 


3.1 x 10 4 



TABLE I: Characteristic length scales for the interaction of ultracold 
molecules in an optical lattice, with values given for RbCs in param- 
eter regimes used in the present work. For the van der Waals coeffi- 
cient, we use C 6 = \M\29E h a\ [20], where E h = 4.3597 x 10" 18 J 
is the Hartree energy and a = 52.918 pm is the Bohr radius. We 
give the confinement length for optical lattice sites with frequen- 
cies (jl>I2ji = 1 kHz and 50 kHz. The dipole length is given for 
RbCs molecules with a dipole moment of \i = 1.0 D, where D = 
0.39343^ao = 3.336 x 10" 30 Cm is the Debye and e is the charge of 
an electron. Here, m r is the reduced mass. 



tic length of the confinement. With the interactions of polar 
molecules having an even longer range, it is reasonable to an- 
ticipate easy creation of a CIR. 

In this paper we study the states and control possibilities 
of two polar molecules isolated in an optical lattice site, fo- 
cussing on the specific example of RbCs. We examine in 
detail the effects of tuning the lattice parameters and elec- 
tric field. We show that the optical lattice can be used to in- 
crease the resonance width past what is possible in free space 
or 2D geometries. We compare the eigenenergies obtained 
for a quasi-2D lattice site to scattering calculations for a sys- 
tem with confinement in only one direction, which accurately 
reproduce the avoided crossings and show their utility for res- 
onant control. Our studies show that tuning the confinement 
has a significant effect on the collisional and bound state prop- 
erties of the pair of molecules, allowing the creation of useful 
CIRs. 



We consider two ground state 87 Rb 133 Cs molecules in a 
cylindrically symmetric optical lattice site. In Table I we 
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list length scales relevant to ultracold molecular collisions 
and give representative values for the parameter regimes used 
in this work. Typical van der Waals coefficients for polar 
molecules are of order 10 5 E^c^ - 10 7 E^a^, much larger than 
those for pairs of alkali atoms (10 3 Ehd^ - 10 4 Ehcfy. How- 
ever, the mean scattering length scales as a oc C^ 4 , giving a 
similar characteristic length to the van der Waals part of the 
potential. We take the dipole moment \i - <fi) z to be the ex- 
pectation value of the electric dipole operator jl for the molec- 
ular ground state in the electric field direction. We calculate 
this electric-field dependent quantity according to the method 
of Ref . [40] . The dipole length, tunable with an electric field, 
is typically the largest length scale in the problem. For a trap- 
ping frequency oj/27i = 50 kHz, = £\^ for [i = 0.16 D. We 
note that the use of a strong dipole moment takes us beyond 
the region of validity of pseudopotential approaches such as 
those ofRefs. [41,42]. 

The molecules are assumed to be rigid rotors, aligned in the 
axial direction by an applied electric field. We approximate 
the lattice site with a harmonic trap and consider only the rel- 
ative motion of the molecules. The combined interaction and 
trapping potential is given by 

V(p z ) = ^(i- —\ + Ql-^L + — — 
r 3 \ r 2 ) r 12 r 6 2m Y p 2 

+ im r (^> 2 +^V). (1) 

Here, p and z are the relative radial and axial coordinates, 
respectively. Also, a> PrZ = 2nf PrZ , where f PrZ are the corre- 
sponding trapping frequencies. The intermolecular separation 
is given by r = ^p 2 + z 2 , and the projection of the relative 
motion along the axis of symmetry is given by m. We impose 
a repulsive short-range potential, Cn/r 12 , setting the Cyi co- 
efficient such that the potential Cn/r 12 - C^/r 6 contains six 
bound states and gives a scattering length of 100<2o- We ne- 
glect the anisotropic C$ coefficient. While arbitrary, setting 
the short range part of the potential in this way allows us to 
conveniently study the important long range effects. 

Although the collisions under consideration involve four 
atoms, the approach described above is justified by the sep- 
aration in energy scale between the chemical bonds within the 
ground state dimers (~THz) and the collision energy or bond 
between them (< MHz). We also note that the van der Waals 
coefficient between polar molecules has contributions from 
the rotation of the molecules as well as the induced dipole 
moments of the electron clouds. An electric field polarizes 
the molecules and changes the rotational contribution. We 
have calculated the extent of this change and checked that it 
does not noticeably change the results presented here, as was 
the case in Ref. [22]. Consequently, we neglect this effect. 
Because we confine ourselves to a single collision channel, 
the resonances we find correspond to shape resonances [29], 
in which the potential experienced by a colliding pair sup- 
ports a near-degenerate quasi-bound state. We note that Fes- 
hbach resonances, in which a colliding pair is coupled to a 
near-degenerate bound state of a different spin configuration, 
are possible for the general case of coupling between states of 
different molecular spin and rotational quantum number. 



We first study the two-body energy spectrum. We solve for 
eigenstates and eigenvalues of the Hamiltonian with the po- 
tential of Eq. (1) using PHAML Version 1.8.0 [43, 44], a par- 
allel two-dimensional finite element code for elliptic bound- 
ary value and eigenvalue problems. PHAML features adap- 
tive grid refinement of the discretized spatial coordinates to 
concentrate the grid in areas where the wave function varies 
rapidly, and high order elements to obtain an accurate solu- 
tion. For these computations we used eighth degree elements. 
Within PHAML, ARPACK [45] was used to solve the discrete 
eigenvalue problem, using the shift- and-invert spectral trans- 
formation to compute interior eigenvalues, and MUMPS [46] 
to solve the resulting linear system of equations. The par- 
allel computations were performed on two nodes of a Linux 
cluster. A particular advantage of the two-dimensional solver 
is its ability to readily account for the anisotropic interac- 
tion and trapping potential. By contrast, an expansion in 
spherical harmonics or non-interacting trap states will strug- 
gle to accurately resolve the wavefunction without a very 
large basis set. However, for analysis of the wavefunctions 
we calculate projections onto these functions. We solve for 
the function F(p, z), where the full wavefunction is given by 
if/(p,z,(p) = F(p,z)e lm< ^. Our bound state calculations con- 
sider only m = 0, but in the scattering calculations described 
below we will consider the effects of nonzero m. We study 
spherically symmetric (co z = co p ) and quasi-2D (oj z » co p ) 
geometries, with the dipoles always aligned along the z-axis. 

The eigenenergies of two RbCs molecules in a spherically 
symmetric lattice site with f z = f p = 25 kHz are shown as 
a function of dipole moment in Fig. la. We use the term 
bound states to refer to those that are bound when the trap 
is adiabatically turned off. States close to the non-interacting 
trap level energies, (n z + 1/2)Hoj z + (2n p + \m\ + l)hoj p , are 
called trap states. Here, Bose symmetry allows the trap state 
quantum numbers to be n z = 0, 2, 4 . . ., n p = 0,1,2,... and 
m = 0, ±2, ±4, .... At zero dipole moment, the trap state ener- 
gies are affected by the van der Waals interactions. Dipole mo- 
ments above approximately 0.1 D cause a significant change 
to these energies. A large number of avoided crossings occur 
as trap states are brought into the potential by the increasing 
dipolar attraction. All crossings are avoided, although some 
are too narrow to be visible on the scale shown. States of dif- 
ferent £ are mixed by the dipolar interactions, with the broad- 
est crossings occuring between states of low I. For exam- 
ple, the first two trap levels converted to bound states as \i 
is increased are primarily of mixed s- and d- wave symmetry, 
whereas the steeply descending state with a series of narrow 
crossings near \i - 0.31 D is concentrated in I - 10. Figures 
lb and lc illustrate the wavefunction near E/h = -50 kHz for 
dipole moments of \i - 0.17 D and \i - 0.498 D, respectively. 
We plot the function p|0l 2 , to more clearly illustrate both short 
and large length scales. For the state at \i - 0. 17 D, vibrational 
nodes of the bound states of the van der Waals potential can 
be seen as rings of constant r < 0.25^ Q . A d-wave component 
at long range provides the head-to-tail configuration, with the 
wavefunction concentrated in the region p < |z|. This makes 
E decrease as \i increases. The state at jj = 0.498 D is much 
more strongly coupled between partial waves, and has a cor- 
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FIG. 1 : (a) Eigenenergies for two RbCs molecules in an optical lat- 
tice site with f z = f p = 25 kHz. Trap states are adiabatically con- 
verted to bound states as the dipolar interaction is increased. Points 
marked 'b' and 'c' correspond to the wavefunctions shown in the 
lower panels, which illustrate the head-to-tail configuration at two 
different dipole moments. For the three lowest energy trap states, the 
partial wave at \i - is indicated. We also give the partial waves 
into which these states have a significant projection after being con- 
verted to bound states as \i is increased. Red crosses indicate the 
non-interacting trap state energies. In (b) and (c), we plot p\if/\ 2 , and 
scale all lengths by the confinement length in the z-direction, l\ o . 



respondingly more complicated configuration. 

We now compare these results to the case of a quasi-2D op- 
tical lattice site, with f z = 50 kHz and f p = 1 kHz. Eigenval- 
ues are plotted as a function of \i in the lower panel of Fig. 2. 
A quasi-continuum of radial trap levels is found, spaced by 
2f p , instead of the strongly mixed states of Fig. 1. Only the 
n z = trap state is within the energy range shown, although 
some bound states have substantial admixture in higher z 
states, as indicated in the lower panel of Fig. 2. The quasi-2D 
configuration has the effect of making the avoided crossings 
narrower than in the 3D lattice case. An intuitive explana- 
tion of this effect is that a sufficiently strong dipole moment 
allows molecules to overcome the confinement which holds 
them side-by-side, and move to the attractive head-to-tail con- 
figuration. The overlap between the asymptotic states repre- 
senting these two configurations is reduced by a higher trap 
aspect ratio. 

It is desirable to attach meaningful quantum numbers to de- 
scribe the states that are observed. This is made difficult at 
large \i by the strong anisotropic interactions; however, some 
approximate quantum numbers can be used. States can be de- 




FIG. 2: Elastic collision rate as a function of \i for a trapping fre- 
quency of f z = 50 kHz and collision energy of k B X 200 nK, show- 
ing the contributions of the m = 0, +2 and +4 partial waves (upper 
panel), and the corresponding m = bound state calculation with 
f z = 50 kHz and f p = 1 kHz (lower panel). The z trap states with a 
substantial admixture are indicated for each bound state not concen- 
trated solely in n z = 0. 



scribed by their projections onto the non-interacting trap lev- 
els with quantum numbers n z and n p , although both van der 
Waals and dipole-dipole interactions mix these levels. We in- 
dicate in Fig. 2 the states which are concentrated in higher 
z trap levels. As we discuss below, it is these states that 
provide the possibility of creating CIRs. States can also be 
described by their vibrational quantum number and projec- 
tions onto spherical harmonics described by the orbital an- 
gular momentum I, which are good quantum numbers in the 
zero-dipole limit. Different £ also become strongly mixed at 
sufficient yu, as discussed above for the states shown in Fig. 1. 

We now make the link between our calculated bound state 
energies and scattering properties in a quasi-2D system. Our 
scattering calculations use the coupled channels technique 
discussed in [22], adapted for elastic boundary conditions at 
short range. We use the potential of Eq. (1) with co p set to 
0, and propagate the scattering matrix in a basis of spheri- 
cal harmonics to a value of r large enough to match onto the 
n z = trap state. The chosen collision energy of &b x 200 nK, 
corresponding to h x 4 kHz, is such that higher z trap states 
are not significantly populated at this separation. Here, &b 
is the Boltzmann constant. We then propagate outwards in 
p, extracting the scattering properties at long range using the 
conventional tools of scattering theory [47]. The results are 
shown in the upper panel of Fig. 2, and agree well with the 
bound state calculations. The resonances at low dipole mo- 
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FIG. 3: Confinement induced resonances created by tuning of f z . 
The top panel shows the elastic collision rate for a quasi-2D trap 
with f z = 50 kHz. We have summed the contributions of partial 
waves from m = to \m\ = 4. Solid lines correspond to \i - 0.3048 D 
and collision energies of x 1 nK, X lOOnK and &b X 200 nK, 
as labelled. The dashed line (ji = 0.3058 D) shows the sensitiv- 
ity of the resonance location to electric field variation. The bottom 
panel shows the eigenenergies of the Schrodinger equation with the 
potential of Eq. (1), using \i - 0.3048 D and trapping frequencies 
of f z = 50 kHz and f p = 1 kHz. The dashed green line shows the 
perturbative result of Eq. (2). Dotted red lines and arrows show the 
intersections of this calculation with scattering states of the given 
kinetic energies. These intersections agree well with the calculated 
resonance locations of the upper panel. 

ment are widest and most isolated from other features, making 
them the most useful for resonant control. Contributions to the 
elastic collision rate coefficient from collisions with higher m 
make the resonance minima nonzero. 

Optical lattices have been used to control scattering lengths 
in neutral gases [38, 39] by making a state corresponding to 
an excited trap level near degenerate with the colliding atoms. 
This depends on the scattering length being of the same order 
as the characteristic length scale of the confinement. With the 
large dipole length characteristic of interactions between polar 
molecules, it is reasonable to expect that similar confinement 
induced effects should occur. In the lower panel of Fig. 3 we 
show this effect for RbCs molecules, calculating eigenener- 
gies as a function of f z while maintaining f p = 1 kHz. We 
assume [i = 0.3048 D, corresponding to an easily accessible 
electric field of 0.67kV/cm. The figure shows six trap lev- 



els with n z = and a range of n p , and a single bound state 
crossing these levels. As shown in Fig. 2, this bound state has 
substantial admixture in higher z trap states. The energy of 
the bound state therefore increases with f z faster than the trap 
levels. 

We have perturbatively calculated the change SE in an 
eigenenergy from changing oj z to co z + Sco z , which results in 

Here (• • • ) indicates calculating the expectation value with re- 
spect to the selected wavef unction. The result for the bound 
state of Fig. 3, evaluated with the numerically obtained wave- 
function at f z = 40 kHz, is shown as a green dashed line. The 
red dotted lines correspond to the energy of a non-interacting 
pair of molecules in the ground trap state, with ^xl nK and 
Jcb x 200 nK of relative kinetic energy. The points at which 
these cross the perturbative calculation correspond well with 
the scattering resonances shown in the upper panel of Fig. 3. 
For a relative kinetic energy of &b x 200 nK, the feature has a 
width with respect to f z variation of approximately 5 kHz, al- 
though accurate electric field control will be necessary due to 
the resonance location being strongly dependent on the dipole 
moment. This is shown by the dashed line in the top panel, for 
a dipole moment of 0.3058 D, which corresponds to a change 
in the applied electric field of approximately 2.5 V/cm. The 
temperature is also of significance, as shown by the calcu- 
lations for collision energies of &b X lOOnK and fexl nK, 
which produce narrower peaks at lower trapping frequencies. 
This result illustrates that the location and properties of the 
resonances can be controlled by manipulating both the elec- 
tric field and the confinement. 

In conclusion, we have studied the role that an optical lat- 
tice can play in controlling the collisional properties of non- 
reactive polar molecules. We have shown that tight confine- 
ment allows for much broader avoided crossings, giving a 
greater resonance width than is available in free space. We 
have also shown that confinement induced resonances can be 
easily created, with the caveat that their location is sensitive 
to the dipole moment. Measurements of resonance locations 
would constrain the short range potential, for which we stud- 
ied just one example with a scattering length of 100<2o- How- 
ever, this should not significantly alter the density of states or 
our finding that RbCs will have several accessible resonances 
for dipole moments less than 0.5 D and trapping frequencies 
on the order of tens of kHz. These results will therefore be 
of significance for upcoming experiments using non-reactive 
polar molecules. 
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